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^ ' Abstract 

We compare flavour singlet and non-singlet vector mesons from first 
principles using lattice QCD. With Nf = 2 flavours of light quark, this 
addresses the LV-p mass difference. Using maximally twisted-mass lattice 
QCD, we are able for the first time to determine this mass difference pre- 
cisely and we compare with experiment. We also discuss uj-p mixing effects 
arising within QCD through the u-d quark mass difference. 



1 Introduction 



The vector mesons are well described by a quark model since they are approx- 
imately 'ideally' mixed, with the flavour non-singlet (p) and flavour singlet (cu) 
degenerate. Experimentally [1] the uj-p mass splitting is small (7 MeV). This is in 
contrast to the pseudoscalar mesons where the rj and rj' are much heavier than the 
vr. From flrst principles in QCD, this splitting of flavour singlet and non-singlet 
meson occurs because of contributions from disconnected quark diagrams. These 
disconnected contributions can be evaluated explicitly using lattice QCD. Early 
lattice results showed that the disconnected diagrams are relatively small for vec- 
tor mesons [2] but are relatively large for pseudoscalar (and scalar) mesons. Here 
we use more precise methods to determine the size of the disconnected diagrams 
for vector mesons and we discuss the resulting phenomenology. 

We describe the vector mesons in a quark model basis for the case of degen- 
erate u and d quarks. In the flavour singlet sector, we then have contributions to 
the mass matrix with quark model content [uu + dd) / \/2 and ss (which we label 
as nn and ss respectively): 

Here m corresponds to the mass of the flavour non-singlet eigenstate and is the 
contribution to the mass coming from connected fermion diagrams while x corre- 
sponds to the contribution from disconnected fermion diagrams. Thus m„„ is the 
p mass. Because of mixing, as will be discussed, rriss does not correspond exactly 
to any speciflc meson. Since the K* meson has no disconnected contribution, one 
approximation is to use the assumption that the connected contribution to the 
meson mass is linear in the underlying quark masses. Then m^s = 2m„s — m„„, 
that is 2mK* — mp, leading to russ = 1.012 GeV. 

The mixing between the nn and ss flavour singlet channels must produce the 
experimental a;(782) and 0(1020). We can express the mass eigenstates as 

mnn + '2Xnn - S; TUss + Xgs + S (2) 

where 6 = Ix^^/ {mss + Xss — fnnn — '^Xnn) to a good approximation for the relevant 
parameter values. The 6 term arising from mixing will prove to be small (around 
1 MeV mass shift) so that one can estimate the u-p mass difference as 2xnn 
and the upward shift of the </> mass due to disconnected diagrams as Xgg- Even 
though the u - (p mixing induces rather small mass shifts, the amplitude mixing 
can have signiflcant effects. For instance the physical (j) meson will have a relative 
amplitude of nn quarks (compared to ss; this ratio is tan 5 where S is the u — (p 
mixing angle) of \/2xns/ {f^ss + Xss — rrinn — 2x„n) which will be of order a few 
percent and can have an important consequence in decays to non-strange final 
states and u production from strange initial states. 
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The contributions x from disconnected diagrams will depend on the quark 
mass. For pseudoscalar mesons, the x values are rather constant at small quark 
masses [3] and decrease slowly with increasing quark mass [1]. Something similar 
would be expected for the vector mesons. One simple phenomenological ansatz 
would be to assume all x- values were the same. Then the uj-p mass difference 
{2xnn = 7.2 MeV experimentally) should be approximately twice the mass 
shift {xss) although with the above model for m<.s this is 7 MeV which does not 
agree. Such quantitative comparisons are not to be trusted for several reasons: 
(i) the p meson is so wide (circa 150 MeV) that the impact of its decay on its 
mass value must introduce an uncertainty of a few MeV at least (ii) the model 
to estimate the mass shift of the 'connected' (meson mass linear in valence 
quark content) is also imprecise (for example, a meson mass-squared linear in 
valence quark content gives 5 MeV instead) and (iii) the underlying assumption 
that the disconnected contributions {x) are independent of mass may be at fault. 
To clarify these ideas, a direct evaluation from QCD is necessary. 

In this paper, we study lattice QCD in the unitary sector with A'^^ = 2 de- 
generate light quarks. This enables us to extract estimates for the u-p mass 
difference. We can also explore the (sea) quark mass dependence of Xnn which 
can help to understand the value of Xss needed to determine the disconnected 
contribution to the meson. 

There is considerable interest in building reliable models of the quark mass 
dependence of the p meson. We show that such models inevitably have conse- 
quences for the uo — p mass difference. We compare these predictions with our 
lattice results. 

We also explore effects arising in QCD from the u-d quark mass difference. 
This causes a violation of isospin and a mixing between u and p mesons. We are 
able to evaluate this mixing (given the quark mass difference as input) from the 
lattice and we compare with experiment. 

Here we follow the lattice methods used in our study of the flavour singlet 
pseudoscalar mesons [5]. The disconnected contributions are evaluated using 
stochastic methods with variance reduction while the connected contributions 
are evaluated using stochastic time-plane sources with the 'one end trick'. 

Use of the remarkable precision obtainable 0, El E] in evaluating disconnected 
contributions in maximally-twisted lattice QCD will enable us to obtain a statis- 
tically significant signal for the vector mass splitting. 

2 Twisted mass lattice QCD and neutral parti- 
cles 

In quenched or partially-quenched lattice QCD, the disconnected contribution to 
the flavour singlet meson does not combine properly with the connected contri- 
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Ensemble 




X T 


/5 


ajjLq 


1% 


ro/a 


Bi 




X 48 


3.9 


0.0040 


0.160856 


5.22(2) 


B2 








0.0064 






B3 








0.0085 






Be 


32^ 


X 64 


3.9 


0.0040 


0.160856 




Ci 
C2 


32^ 


X 64 


4.05 


0.003 
0.006 


0.157010 


6.61(3) 



Table 1: Summary of ensembles produced by the ETM collaboration used in this work. We 
give the lattice volume x T and the values of the inverse coupling f3, the twisted mass 
parameter a/iq, the hopping parameter k = (2amQ + 8)^^ and the Sommer parameter r^/a in 
the chiral limit from ref. [TTl[TUj. The data sets cover 5000 equilibrated trajectories (10000 for 

bution to give a physical state. To avoid this problem, it is mandatory to study 
flavour singlet states in full QCD - with sea quarks having the same properties 
as valence quarks. Then the spectrum of flavour singlet states is well defined 
and can be extracted from the t-dependence of the full correlator. Here we focus 
on the case where there are two degenerate light quarks (called Nf = 2) which 
is a consistent theory in which to study the flavour singlet mesons. We use the 
twisted mass lattice formalism [3, E], for a recent review, see ref. |9]. 

The results presented in this paper are based on gauge configurations as pro- 
duced by the European Twisted Mass collaboration (ETMC). The details of the 
ensembles are described in ref. [10] . In particular, we concentrated for this paper 
on the ensembles labelled Bi, B2, B^, Bq and Ci and C2. We have compiled the 
details for those ensembles in table [H The S-ensembles correspond to a value of 
the lattice spacing of about a ~ 0.09 fm and the C-ensembles to a ~ 0.07 fm. The 
spatial lattice size is of about L ~ 2.2 fm for all ensembles used here, apart from 
Bq, which has identical parameters to Bi but L 2.7 fm. The corresponding 
pseudoscalar mass is in the range 300-440 MeV. 

The twisted mass formalism at maximal twist is order a improved and allows 
access to light pseudoscalar mesons [11]. There is one complication, however, 
namely that the twisted mass lattice formalism breaks flavour and parity sym- 
metries at finite values of the lattice spacing a. These symmetries are restored 
in the continuum limit and the theory is well defined (so providing a valid reg- 
ularisation) at finite lattice spacing. One consequence of this flavour-breaking is 
that charged and neutral non-singlet mesons can have a mass splitting (of order 
a^). Moreover, the neutral non-singlet mesons have (order a^) contributions from 
disconnected diagrams [121 [H] . Note that the order splitting between p"*" and 
p° was found to be compatible with zero [3]. 

In the case of twisted mass fermions, the bilinear operators appropriate to 
create the u state are OY = ip'yiip and Of = ip'yi'y^ip which on transformation into 
the twisted basis (used in lattice evaluation) will become xiiX and X7i7475''"3X) 
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respectively, where acts in the (m, d) flavour space. This latter case amounts 
to evaluating, in the lattice basis, the difference of the disconnected loop between 
u and d quarks. This enables a very efficient variance reduction technique [HI El [5] 
to be used to evaluate the relevant disconnected diagram for the case of the 
operator. Here we follow in detail the methods to evaluate the disconnected (D) 
and connected (C) meson correlators described in ref [S]. Some discussion of the 
connected correlators for vector mesons has already appeared p!H| [T3] and a fuller 
study is in preparation [T3] . 

The connected meson correlator describes the flavour non-singlet meson (p), 
while the combined correlator describes the flavour singlet (u). Assuming suffi- 
ciently large t-separation so that the ground states dominate, then we have (for 
Nf = 2) 

= cexp(-m(p)t); C{t) - 2D{t) = dexp{-m{u)t) (3) 

Hence 

2D{t)/C{t) = exp(-At) + 1 (4) 
where A = m(cij) — m{p). Moreover, if A is small, as we shall find, then 

2D{t)/C{t) = ^At + {l-^) (5) 

Since the difference between p and u is so small, one would also expect the 
creation amplitudes to be equal so that d c. Then keeping the first order in 
the mass splitting, this leads to the simple result that 2D{t)/C{t) = At + 0{A) 
which will be a good approximation provided At -C 1. 

In twisted mass QCD, there is an order contribution to the p° meson prop- 
agation arising from disconnected diagrams. This has been studied previously [3] 
using the operator which enables the variance reduction method to give pre- 
cise results in this case. No statistically significant contribution was found. For 
our present purposes, however, we wish to focus on the disconnected contribu- 
tion to the u meson compared to its connected contribution which is in turn 
equal to the connected contribution of the p° propagation. So we compare the u 
propagation with that of the connected contribution (C) for the p° meson. 

We present some results for this ratio 2D /Cm fig. [TJ We have evaluated this 
for local meson operators (L) and non-local (fuzzed) meson sources (F). 

In order to interpret this signal we evaluate the effective mass from the con- 
nected component which is displayed in fig [2j This shows that the ground state 
dominates in this case for t/a > 5 for the LF and FF cases. This allows us to 
interpret the data in fig. [T] as being dominated by the ground state for that t- 
region, so that one can determine the mass splitting from the slope as described 
above. 

We fit the behaviour of 2D/C as a straight line (eq. [5j) in the range of t/a 
from 5 to 12 (6 to 12 for the LL results) for (3 = 3.9 and 6 to 16 (7 to 15 for 
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Figure 1: Ratio of disconnected to connected contribution to neutral vector meson (created 
by O"^) versus t. Results are from ensemble Bi [L — 24) with ajiq = 0.004. The ratio is given 
for cases with local operator at source and sink (LL) and for cases when one or both operator 
is non-local (fuzzed F). The line represents the fit to the LF data in the t/a-range 5 to 12. 



LL) for f] = 4.05. The of these fits is acceptable. We take the average of the 
slope from the LF and FF fits and we determine an error for the slope aA from 
combining the range of values found (including the LL case). The intercepts are 
very close to zero, so the impact of the factor d/c on the determination of A in 
eq. [5] is negligible. We find that the results from two different volumes (spatial 
extent L=24 and 32 at fiq = 0.004) are consistent with each other. 

In order to compare the results from different lattice spacings, we plot them 
using the ro/a- values given in table [T] to convert our results to GeV units using 
ro = 0.454(7)fm obtained [H] from the ETMC evaluation of J't^. We show the 
uj-p mass splitting in fig. [3] and we see that our results at different lattice spacings 
are consistent with each other. 

Extrapolation to the physical limit (for u and d quarks) gives m{uj) — m(p) = 
27(5) MeV assuming a linear dependence. Since there is limited evidence for 
the form of the extrapolation, we assign a error which covers the possibility of a 
constant extrapolation (namely 27(10) MeV). 

This may be compared with the experimental splitting [1] of 7.2 MeV. In our 
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Figure 2: The effective mass (in lattice units) of tlie connected contribution to tfie neutral 
vector meson (created by O^) versus t. Results are from ensemble Bq {L — 32) with a/ig = 0.004 
and rg/a — 5.22. The result is given for cases with local operator at source and sink (LL) and 
for cases when one or both operator is non-local (fuzzed F). 



QCD study, we do not include electromagnetic effects or effects arising from the 
u-d quark mass difference. These effects are estimated [16] to be of the order of 
a few MeV. Furthermore, the definition of the p mass is uncertain to a similar 
extent, because it is such a wide resonance. Hence precise agreement between 
the observed u-p mass difference and that found from lattice QCD is not to be 
expected. 

The quark mass dependence of vector meson masses is of considerable inter- 
est, especially in the lattice community. Next we discus models which may be 
compared with our results. 

3 Models for vector meson masses 

The disconnected quark diagram that contributes to the u — p mass difference can 
be used to inform models of the mass- dependence of the vector meson masses on 
the underlying quark mass. Such models contain contributions from two-particle 
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Figure 3: The uj — p mass difference versus the quark mass given by mp. A linear extrapolation 
is plotted. The lattice spacing is a and the spatial extent of the lattice is L. 




Figure 4: Disconnected and connected two-particle intermediate states. 



intermediate states evaluated in some effective theory. One contribution to the 
u—p mass difference comes from such two particle contributions to the self energy. 
This gives an excellent way to check on these effective Lagrangian models for the 
p mass versus m^. Here we retain the physical meson names (p, u, it) to describe 
mesons with these quantum numbers as the light quark mass is varied. 

One way to visualise this contribution to disconnected diagrams from two- 
particle intermediate states is illustrated by the two quark diagrams describing 
disconnected and connected quark line contributions in fig. HI 

The two-body contributions (having non-analytic behaviour with m^r) that are 
usually considered are the pseudoscalar-pseudoscalar (PP) and vector-pseudoscalar 
(VP) intermediate states. Here we tabulate the relative contributions to p and u 
mesons: 
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V PP VP 

p\ IT TT UJ-K {Ig'^) 

u: - pTT {3g^) 

The notation for the VP case indicates that — > pvr has 3 terms {p^ir'^, p^vr^, 
p^7r+) whereas p — > ojtt has just one {uj^tt^ but with the same couphng). Thus 
the contributions from these intermediate states to the uj — p mass difference are 
VP{2g'^) — PP while the p mass has VP{lg'^) + PP. This exemphfies the hnk 
between models for the quark mass dependence of the p meson and models for 
the a; — p mass difference. 

We now discuss some models proposed and focus on their implications for the 
uj—p mass difference. We first consider the VP intermediate state. The behaviour 
of this contribution versus quark mass is as and estimates of the contribution 
(without using any form factor) to the p mass at the physical pion mass are 
3.4 MeV [in] and 5 MeV [TTj. The strong dependence on implies that, in 
these models, the VP contribution for a pion mass around 300 MeV will be much 
bigger: for the uj — p mass difference they obtain values of -70 MeV [I6] and -30 
MeV [TTj using a form factor cutoff. This large value obtained without cutoff 
is quoted as 'not to be taken too seriously' by ref. (TB]. For more discussion of 
regulators which may improve the convergence, see also ref. [IH]. In our opinion, 
the evaluation of the VP intermediate state is delicate since as the quark mass 
decreases, the V and VP states become nearly degenerate. In this circumstance, 
one needs to take account of the width of the p meson which could modify the 
expressions significantly. 

The PP intermediate state (i.e. vrvr contribution to the p) has a non-analytic 
component behaving as m^log(m7r). One approach [l7] uses the relativistic ma- 
trix element but introduces a form factor to emphasize the kinematic region where 
the model should be more reliable. Fitting to lattice data, they obtain a contri- 
bution to the p mass of -35 MeV for the physical pion mass and -40 MeV for pion 
masses around 300 MeV. These values imply a contribution to the uj — p mass 
difference of 40 MeV for pion masses around 300 MeV. 

In summary, models give quite large contributions to the uj—p mass difference, 
with the VP intermediate state dominating, and becoming more negative as the 
quark mass is increased. However the models have an increasingly large contri- 
bution at larger masses which needs to be tamed by a cutoff or other method. 
Our lattice data show a decrease with increasing quark mass, but with numerical 
values that are much smaller than those typically given by models. 

4 UJ — p mixing 

Within QCD when m„ ^ rrid, isospin is not a good quantum number and there 
is a mixing matrix element between the a; and p mesons. One experimental 
manifestation of this is the decay uj tttt which creates an interference pattern 
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in the dominant decay p — > tttt near the uj mass. The mixing can be defined from 
a mass matrix (with basis states {u'-)iU ± d'yid)/\/2 ): 



(6) 



Phenomenologically [I9] the transition matrix element T^p = —3.1(3) MeV. Some 
of the observed effect can come from electromagnetic contributions and these are 
estimated to contribute between 0.4 and 1.2 MeV. So the QCD contribution 
would be around -4 MeV. 

We can study this contribution on a lattice by considering the cross correlator: 
create a cu meson (using (ujiU + d'jid) / v^) and destroy a p meson (using {u'jiU — 
d'-iid)/^/2 ). This cross correlator will have a connected piece (like (C„ — Cd)/2) 
and a disconnected piece. The disconnected contribution simplifies using the 
identity D^d = D^u yielding {Duu — Ddd)!'^- Here D^d means the correlator 
between a w-quark loop at source and a d-quark loop at sink. Because of the form 
of the cross-correlator, it may be expressed as a difference between propagation of 
vector mesons containing quarks of different mass. It may then be related (using 
techniques similar to those used in ref. [20] and assuming that (m^^ — mp)t <^ 1 
to simplify the expression) to the mass- matrix transition element: 

_ (krip^ dm^ m„ - 

where we have assumed that we may retain the first power of m„ — only. 
Here the derivatives are with respect to the degenerate quark mass in an Nj = 2 
theory. These derivatives are those that we have access to from our study above. 

We can estimate this transition rate phenomenologically, using drrip/dm'^ f» 
l/(2mp) (from the (j), K* and p masses, see ref. ^T\) and dm^/dniq = m^/rhq 
(from lowest order ChPT). Here rriq is the average light quark mass (m„ + md)/2. 
This gives us 

and with mu/md = 0.56 from lowest order of ChPT [1] we obtain T^p = —3.6 
MeV which is close to the experimentally determined value. 

We now discuss evaluation of the above expression in eq. [7] using lattice results 
from first principles. This involves the quark mass dependence of the p and a; 
masses in an Nf = 2 theory, which we have available. In particular we are able 
to address the quark mass dependence of the disconnected contribution to the uj 
mass, which has not been studied previously. 

Since drrip/dmg + dm^^/dmq = 2dmp/ dntq + (i(m^ — mp)/dmq we compare 
these two terms. We can evaluate (i(m^ — mp)/dm?p as the slope shown in fig. [31 
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We obtain a slope — 0.15(15)/(2my) where the error encompasses the possibihty 
of a constant behaviour which is not strongly excluded by our lattice results. 
Compared to the phenomenological estimate that dmp/ drn^p l/(2my), this 
value indicates that the contribution of the uj meson disconnected diagram to the 
uj — p mixing is relatively insignificant (at around 7% )- in agreement with earlier 
phenomenological estimates |16j . 

We find that lattice QCD naturally produces effects of the correct size to 
explain the observed uj-p mixing. 

5 Summary 

Previous lattice results were very exploratory: for Nf = 2 sea quarks of mass 
(given by {rompsY = 3.7) corresponding to about 1.5 times the strange quark 
mass, an uj-p mass difference of 2(3) MeV was obtained [2] with quite strong 
evidence that the sign of the effect should be positive (the same as found here). 
A more recent study finds [22], in the chiral limit, an u mass of 790(194) MeV, 
fixing the scale at the p mass: so a mass difference of 15(194) MeV. 

Here we have presented results coming from a well established signal and we 
are able to address issues such as extrapolation to the continuum limit. We see 
in fig. [3] that our results are consistent with each other at different lattice spacing 
and for different lattice spatial volumes. There is evidence for a reduction in the 
splitting with increasing quark mass. Since there is limited evidence for the form 
of the extrapolation to the physical point, we assign a error which covers the 
possibility of a constant extrapolation (namely 27(10) MeV). 

This may be compared with the experimental splitting [1] of 7.2 MeV. The ex- 
perimental result includes electromagnetic effects (estimated to be of order a few 
MeV) and the definition of the p mass is uncertain to a similar extent, because 
it is such a wide resonance. Hence precise agreement between the observed u-p 
mass difference and that found from lattice QCD is not to be expected. What we 
do establish, however, is that there are substantial contributions to the u-p mass 
difference arising within QCD from the disconnected contributions (known phe- 
nomenologically as OZI violating contributions). Moreover, these contributions 
are of the same sign as the experimentally observed difference. 

Our study shows that it is possible to extract a signal for the u — p mass 
splitting versus quark mass and that this quantity will be very useful in refining 
models for vector mesons near the chiral limit. Indeed some existing models give 
results which are significantly different from those we obtain (especially for the 
VP intermediate state contribution). 

We also showed from first principles that the observed uj-p mixing arises nat- 
urally in QCD with the correct magnitude. We were able to estimate, for the 
first time, the contribution to this mixing arising from disconnected diagram 
contributions and we showed that it is relatively unimportant. 



10 



It has been suggested that the uo — p mixing is the mechanism behind charge 
symmetry breaking in nucleon-nucleon interactions. However there are other 
competing mechanisms (see Miller et al. [19] for a review) that try to explain 
charge symmetry breaking. One problem with the proposed mechanism oi u — p 
mixing for charge symmetry breaking in the nuclear force is that the momentum 
dependence of the mixing was hard to constrain. Although we have only evaluated 
the uj — p mixing at zero momentum, our techniques can in principle be used to 
investigate the momentum dependence [22], so our calculation is the first step 
in determining the mechanism of charge symmetry breaking in the nuclear force 
from first principles. 

Although we have studied QCD with Nf = 2 light degenerate fermions, one 
can extend the discussion to include the strange quark by assuming that the 
disconnected contribution is independent of the quark mass, so yielding an esti- 
mate of Xns ~ Xnn = 13(5) McV. This implies that the (p-u Hamiltonian mixing 
amplitude given by V2xns should be of comparable size. This is in qualitative 
agreement with phenomenology [16] which finds values of 8 to 14 MeV from an 
analysis of (p decays. 

The mixing between the u and (j) is required to understand B and Dg decays 
with an or in the final state [211 EH] • For example Gronau and Rosner [25] 
predict the branching ratio B{Df — > uje^Ve)/B{D'^ 0e"^z/e) in terms of the 
u — (f) mixing angle. Other mechanisms can contribute to this ratio, such as 
"weak annihilation", so a first principles calculation of the u — (p mixing angle is 
valuable. 

In summary, the vector mesons have a rich structure beyond 'ideal mixing' 
and this rich structure can be evaluated accurately using lattice QCD. 
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